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(57) Abstract: An image 
processing system is described 
in which sets of image elements 
having display values outside of 
a target range B of display values 
are each respectively morpholog- 
ically dilated. The intersection 
between the morphologically 
dilated sets of image elements 
is then identified and- those 
image elements within the 
intersecting region are removed 
from the set of image elements 
having the target range B of 
display values. This removes 
image elements incorrectly 
appearing to have display values 
corresponding to the target 
range B of display values due to 
aliasing effect between regions 
of image elements having 
display values either side of 
the target range B of imaging 
may be two-dimensional or 
three-dimensional imaging. 
The morphologically dilatation 
is preferably performed 
with a quasi-circular or a 
quasi-spherical structuring 
element having a radius of 
between two and three voxels. 



JP 2004-525715 A 2004. 8. 26 



(i9)B*@«iTjf<jp) (12)45 ^ 4t Wt & SI(A) (li)wwajffl^s^ 

4fft2004-525715 
(P2004-525715A) 
(43)&SB V*16#8326B(20O4.8.2e) 



jmj irit . trl. 


F 1 








a o i q a /no 
AO 1 D O/Ua 


A6 1 B 


6/03 


3 6 O J 


2G088 


MOID O/UU 


A6 1 B 


6/03 


3SOG 


4C093 


Ao i d 5/055 


AG 1 B 


5/00 


G 


4C096 


AO 1 □ O/UU 


A 6 1 B 


8/00 




4C60 1 


CO 1 R 33/28 


GO 1 T 


1/161 


C 


5BOSO 






(£40K) «BSHC»< 


(2i) 


*fJS2002-582461 (P2002-582461) 


(71) &® A 


500263460 




(86) (22) UJH-B 


¥^14^2^210 (2002.2.21) 






■J H 


(85) SKJCtittJB 


¥m 5^10 200 (2003.10.20) 






(86) bmhi*** 


PCT/GB2002/000768 








(87) 


W02002/084594 






(87) 


¥^14^10^240 (2002.10.24) 




<#tfc£L) 




(31) Cfttt±g»4 


0109524.9 


(74) ttMA 


100064908 




(32)9*0 


^13^4^180 (2001.4. 18) 




#«± 




@3)«5tsflfc£aa 


KB (CB) 


(74)ttJ®A 


100108578 










#JS± »« 








(74) rmx 


100089037 










@& 


■ 






(74) 


100101465 












iEffi 













(54) HIUettR] H«r-^Sa3lC^»5Sfli<DA»Wft«jt©»iE 



(57) imn 

n»4, BI#OSS»BfflBcDi^{RiJcD«^ffl%W-r§iii#S 



Start 
T 



10 



| Capture voxels characterised by density D 
| Identic set of voxels IV with Bn* a;Oiai», [ _X~ 12 
| Identify set of voxels D„ with Q>R TO \-^~^ 



j Identify set of voxels Di. with D < B™» 



Apply approximate spherical morphological 
dilation to Du and Dl to form D'„ v D' L 



K 1. 
18 



Identify set Dp„ given by Intersection 
of O u and I> L 



Subtract from Dut to form Dv 



Display DU 
i - 



-20 



End 



(2) 



JP 2004-525715 A 2004. 8. 26 



[ PF » * <D IB B ] 
C»#3l 1 ] 

K is ^ T „ 

B m , n 6 B m a x $T?©*wffi©S«lSHrt©**ffi*W*SB«©l!fltM©IH* 

k m * S x t- v 7 t . 

M!E»lr->«jSffi©BttgSK©ffl»C» J&«±©»3B*»fl! LT. fig® L ffi ^ ffi © B4R 

m m s m © m m. m m © m a* , mm&t> o (ommmmnmoytpic zwmw&*%i 
c m * m 2 ] 

§IS*m l tEl«SSo 

c Sit** 2 fciBa£©73 2£o 

4 ] 

imxm 5 ] 

3 $ fc » 4 tc 12 IS © 7? r£ o 
[ It # « 6 3 

is # * 5 k re « © # s o 

[ » # * 7 ] 

3ftSM3afi©S9£ffi**LT^Sc t^^fttt-^W*^ l 6 6C5^0^fn* 
- « tc IB « © 73 r£ „ 

[ m * m 8 ] 

itu f B ft # & S © $J £ {I f i , + MRIX^t>, i^SX + t>, PET© 

3^©io*ffl^T«fflSft*cfc*#«fcT5»#«7fcJE*©73i£. 
[ IS * ^ 9 3 

MIEHB^flOlIifcgJROffltiU 31 p> (/•> iB SI {c ffl ^ L , iff IB ilf {fi © SHI* g 3g © ffl f£ , # 
K: ffi S L , ft»-o[t(lIS@«©H«EJR©llfl{4. 3>'h7XbtfM;&-&A,T^3ifo<§k:tgS 

tic t««F«i:-rsai#5i i *» e. 8o??ot>fnfr-«KEi©ss, 

[ W * 91 1 0 ] 

S^A^^fI^r*LTi/^iafiS«©E^Jlc < l:-pT^fig^tifcHft7 :? -^^:®S-r?)^B 
*5 <^ T , 

B m , „ U T © & ^ ft * * T § fg i> * 55 m © H ft £ * © ti % K 9J T 2> » ft rT &g * ffi ^ S 
^M©ffi©^SUg§i:. 



(3) 



JP 2004-525715 A 2004.8.26 



huei« v^*^ffloB«S3isoiHic, ^ffiioii^ifflUT, m^Lrcs^mmmvmm 

M^O*fc#ftrSSt>DOiB«BJRfl!)ffl*»giJ-rS»ff^RrfiB!a:S5*>»)WgiJ«i:, 
MSEB«<OB«gJR<Dfflfr&. M IE X to D © B « S JR © ffi <0 + tc t» # & ■ « g * * 
D US «^ T , «iE*nfcB«oiii*fi3R<oa*ff»)a-rttf^3b<Brfil*XteOI»*»i:*«il io 
TV* C ^fllttSgio 

1 1 ] 

ffi mm mm m « * ? -a a* t* & d , uEs?jii**-ir;K0£*ssM?*8c t * w ® t -r 

[ » * Jfi 1 2 ] 

^SCfc^^ittSISfSl 1 lc IB ifS co 2§ B 0 

im&m i 3 ] 

t^Ck*|$lif5||*Sl 2ICGIOS1. 20 

1 4 ] 

i 2 $ fc a i 3 ie « co m m . 

C » * Jg 1 5 ] 

m * i 4 k t a * co g e . 

C » * * 1 6 ] 

n - « la m co s a o 30 

[ » * « 1 7 ] 

to f E ft # 3£ g © M fit f± > CATxtt>, mr ix + t>, ifftx + ty, pet<o 
[ m $ 9i i 8 ] 

SttlEffivffi©B«g*©JHH:, S&frt^fiftlcffiliiL. ttilBK WI © H 1* g HI © IB « » # 
tUSL, o m-e g«R®Bft£5Re>ili§tt. a > h ^ x h & &t*^ & me taam 

[ » * « 1 9 ] 

B m , „ B m a x *^©£^tt©l«IEHrt©«jSffl*#-rSB«<DB«Hfli©ffl*' 

^ fit CO ffi CO fi gij D V s y ^ t , 

V * jS fit co ffl co IS 9J D V v 5 t. 

M IE fg v 3! ^ fit © B # g SS © ffl K . JB!»±©B3S*iIfflLT, K9l/ii^SijiiOifll 

M IE « V St ^ {§ © B ft g SI © ffl tc . M±©M£iifflLT, KfiL/iilV^^ffiOitt 50 



(4) 



JP 2004-525715 A 2004.8. 26 



ffiK»SBLfcffiv£iSli©ffl«B*<Dfi:fe <fc M IB M L /ci&i/>S^ffi©ili{i!gfR©;SS© 

im&m 2 o ] 

«} IB H # g IK fi * -tr ;b T* & 9 . tt IB iE 5>J t± # * -fe ;1/ © = & tc E 5U 7? 5 c i: * * 1 1 
[M$«2 1] 

[ Bl * S 2 2 ] 

i/^5C t^^f tt§ll*S2 0 S ft ti 2 l tcK«©nyif a-^ 7Dy5LSn D po 
[ i ft * JS 2 3 ] 

2 1 Iftti2 2{cfB«c<D3>'trj.-^yD^^A^p c So 
C«*®2 4 ] 

mmmWB(DMf$.mmi*. te v> tc ^ 2 . 5#^t*C¥I%SLT^4i:l:«Wti:t« 
»*«2 3 CfBi03ytfa-^yn^5Aiip a p 0 
CW*S2 5 ] 

^nSffltaSEOSISf *SbT^4i:i:*1tai:t5B«ai 9 £2 4 © ? © <^ "T 
[ W * W 2 6 ] 

WJ IB fit 3S K <D «) j£ ffi fck , C A T X + -V > . M R I X + ^ V , jS ^ X > > PET© 

3 "5 © 1 Ofcffl^T&tBStxSC fc^r#^fc-r?,S»^«2 5 (CIS* © 3 > M zl - % - rf P 

2 7 ] 

lWKffi^ffi©H«BfR<Z>ffltt* £ £ ^ « fc ffi S L , W IB SI vW« © H (ft g * © £ It , # 
tStt^^i.ttsBf^i 9 ^ p> 2 6©no^fn^-lt:ffil©3yifa-?- 

[ fg OB © & US ft IK BJ! ] 

[ 0 0 0 1 ] 

* #g m it , Hflii^-^ffiao^wteH-rs. 1$ tc , # ?g bj§ t* , s^ffl©<fc9j£v>:£®Brt 

lctiLWt%. S^i©fe5IBlrtfcS5SSI%*t§I«SI ( image element) © * 
^ IC ffl T & o 

[ 0 0 0 2 ] 

— &7tmm*&? t><DT2bti, =*7cH«*^-rt.©-e*n, h & g sn # , 

fiMSfcJf9-««^ffiK«i;T, ffi6nftlig|©'>4< i: & < o © ill iR *5 <£ t>* 0 IE 
„ fe?)«-g-{C«, £ tlP><DM&&Mlt. ^Pfi©&-g>^H (defection device) S> © {I -!§ 

si; ft © $y £ ^ a l t , is, ft#©i&s«fc«rtetei»3!ttw-&na. - 0»j t l t . e ^ 
fflai^'NwiSffl tc*3 ^ t, c a T^ + t ysftSMR i x + ty^6i?n«f§ 



(5) 



JP 2004-525715 A 2004. 8. 26 



3/cftfC, I^Sfe^S/^f <0S*5SitCYy tTV^LT, #5i£©f#®, fflxlfif 
[0 0 0 3 ] 

otSJti5> 3>h^xhif33iSij3:-&A,T*^£>ifti<g6#, c©B#rt{c#&t-&o 

[ 0 0 0 4 ] 

T^S. & ft , «©«tt2fc*&a*vfi»4£©IWK:, #©SRlg2©igiV>ffii:*£fr 
5 o ^ < © *1 ^ . ctx6©^^©iiilR(i. rtie6K:ffiMi*-S3SjSffli:a«£Hi;*^ffi#*r 

ltv^ 0 tti«>^.te, I2oii^, »i«aBt«ft»ic, 2 *> & it ft fa e> © # 
s t? jfli s 6 t ffl s -r a « ^ « * w -r * a ^ s * * a «: ft k «i a -r * t , istss 

Sfi§ft5J:?Ki^5H' l J7J/>'^ (aliasing) illCtoT, g? o T , #©«g^ 20 
[ fg W © P*fl * ] 
[ 0 0 0 5 ] 

B m i n 3&» 6 B m a , *-P©a3Sffi©B««Hrt©*^ffi*#-rSB«©H«BiR©li* 

B m a x & v & m^£^m<Dmmmm(oM*m%\-f & xt- v t , 30 

M IE {g {i © H ft! S #ft © ffl te , JB«B±©»3B*afl!LT\ illfcIi/^ieoHi 

MIHKv>a^fiI©iH#^m©)H»c, JB»±©»31*afl!LT, III fciSv^^ fit ©B& 

mmmmi,rci&^mmmcDmmmm<Dm^ ^xsmmm^L rcm^m.mm<Dmmmm(Dm<D 

p^©4'tc#ft-rs, 2fcftt>©!B«H*©ffl*»8'J-r3;*-rv:/i:, 

WEB«©iii«S*©ffl^6» MK£fcD©fflffcg*©IS©*fcfcSF&-rSHi«gilS«:8l 
[ 0 0 0 6 ] 

# fg Wtc T « , Iffl » o T\ S^l<DlISSi(cffiit§SiI^tt§ <t3tclx.§ 40 

, f n6*ffiDH!< C ttf-e*5. gM1K«©iiiofl!lfc*Jt*S«« ( -f & ft a 

«i S Ift © ffl ) *«»WK»a**T, J tft5©BB»L;fcffi«©3Efc9*#ftS£:i:fc:<J:o 

tXbiOT, C ft fS i; T , COfitdtt, ±fB©«H?««*aJRW»cK»JLT, BSbS 
iS?i:g^fi5H^©«^{iI^#-raiii^S«©««t©P^©«^tcfeaiii«S*^J&oTL 
cft5©*#©IH«E*# t Ht8'J£ftfttf, ^ ft e> £ , g«RIGHF*i©£CT 

ftfcSH*RjR©ffifc, ^W©x-r .'J 7 -> :/^.© A&tt£«ii (artefact) ^IS^ttS 50 



(6) 



JP 2004-525715 A 2004. 8. 26 



[ 0 0 0 7 ] 

# fg B£ © & ffi fi , iEL<KB«J*nfcB«H*©**K:a^:«:IB«*#iiaci:*<, o 

4«*fcttaia5fc:t,±»*L*3i»*s. m * , n 3b« ^ % hb * a . h 

[ 0 0 0 8 ] 

£7&3? 0 :£ f£ HJ3 & . H^7C#^-b;l/ (voxel) x — 2 <D M <D X m ^ % <D IC, 
[ 0 0 0 9 ] 

mm±<D&m (dilatation) fct x M 4 <D # U Ic fS C T > C ttfT'?5 

*P (quasi) J*«©flH©««F»i©£T©j6fc:3g&Sn<5. 
[0010] 

C comB±(DB^iC m V h ft &1$&<Dmi&mM (structuring element) It, V 3 V 3 & 
A * £ % m ct^TtSA', * % BJ! {c v> T (± , C <D #t l£ K M ft . 2 e> 3 sl< * -tr ;V 

© a ^ £ © ¥ m * * t z> t * «fteja«wi?a&D, *e ^ fc v 2 . 5 # -tr ;b © a # * r* 

SSCkiiC, <fc D £f £ L ^ C £ # # 2r> o tc o 
[001 1 ] 

& x x A tc *5 ^ T ti , * ^ ffl fi , as, i»W*#tti:*5:Jn§t(?), « * €l * 
fi5Sfit{±*ffSLT4JP)-r, ^ © ft *> t> . S!l S X ^ A , l^lfCATX^tt, MR 

o 

[0012] 

* m wit. mmmm. mz.ttjbi'gmmvmmftv, v ^ ^vm&mtt^ t^x^^^m 

mbfr^Mi® t . # £ # E gij T 3 C £*Ufr2>miC WMfrZ A2t> W & #§ ig Bfc* "T 
& © {c , »ICJ;<1LT^5. 
[0013] 

[0014] 
[0015] 

H! « 5* - * , «aU(fCATX + tt« M R I X * + A , gfSX + ttSftliPET (Po 
si tron emission tomography system; P§<! : f#{*t®T/lJJiK> / XAA) 7 1R ^ £ tl fc 1 0 
0 ft © 5 1 2*5 1 20-^7CiB«0*a*» UftOgftK ±5i«fflIC^lJT, J» & 

MCC c T IK B£ L fl: ^ „ 
[0016] 



(7) 



JP 2004-525715 A 2004. 8. 26 



s?$ag©ft46©g^ffi©i&B£-as?-r & c tt>K aetin^. -mht, f$s©®a 

[0017] 

LTiLftUF^^AtJS. @3Kfx L/ciat, f|i»tf;&3g«f©5BB(t, 3. — If 
fc,£-s>TiBS?£ft&±RgB m a x *5 £ If T PS B m | „ a. — if it, i5 ^ 6 < 

ifoeiS«Bi*2lfTrSi:*»cg#SftS«fc5ax =r v h ^ x h*3aay*-&A/TH**ilii** 

{i © J: t)iBJ£v^iiSHrttciB«*n> B«ffiBB©TffiiJ©«?ii£::&-g>T<D©{I^{I©®B 
Afc^oTftOWSft&c H«fC, iS ^tJ©fii©®B C ft, ittffiHB<OffiV{RJ<O0tl?& 
ilfixB©J:?fc, ttBB A rtom^ffi*ig-rffl«^, SB B C rt © ft *§ fit * M * fi tK fc * 

ft § o 

[00 1 8 ] 

or*, m&&&is.-? %>rctb<D#2 -tfrtfi. *t js -r ^m^m^^-r * t&icffiis.-$ ns 20 

^5. Xr'?7l 2 1?, H 3 fc ^ L ft , ^--■9 t ^^^Lfciifi{i<DB^®HBrt<DiiS^r 
St5t^-b/KDlA^9J^n§o X t- -y 7 1 4 T' . B 3 0«lCKiStSl«*f t 
5 * -fe ;U CD ffl «t M BU * n « o X t- y y 1 6T, B3C«fflAli:flatSifi*gt5# 
* -fe ;b © ffi # m m t£ ft S o 
[0019] 

x r- -y 7° l 4 <fc tf l 6 T ! 18 gij £ ft ft jK * -b ;l/ © 3S ft , * ft *"> £ , * ft ^ ft © IS J£ © J£ fS 
± © K 31 ic frit Z ft , »ffiLfcsK*-fe>l/©l|§#£j££ftS. J£ ffi ± © B « ft . ^5^5* 
»*i5Ct*«T?t5*', C©0ljT*ft, ^'\©#^-tr;l/jfiffilJCSt3<^^m© (quasi-sp 
herical) ISSl^ffll^. C © E£ ft , 2 ^ 5 3 * ^ -b ;K0^f 5: t L T ^ S »S 30 
L < ft , ft ^ ft ^ 2 . 5#^-b;l/-e*S. C©£9&J*ffiJ©«J3ScT-©JEJ^±©|jg?I(t, £- 
3R^-b;bffl*, t;^BA,t'l>T, © « j£ fc <fe o T £ « * ft 3 titi* F*3 © £ T 

<D#*-fe;Hc^ta-r*«fc3K:f^ffl-rs. eft ft, H«LT^SW«*t>-r^fcl£tf*/* 

31 $ -a- s „ 

[ 0 0 2 0 ] 

Xf •?7 , 2 0f, 20©^3i^-&e>ftftil^S*©)Sa% Jt'« lt> M75-©H©4>fcgift5 

wwtwmtk^m? & o cfte,©ii#s^{t, 2-D©ffl©r B i©«^^s«{-*@^-r^.o e ft e> 

©8t#««tt, # * *A*iBl5 Tx-f U T S" ft § (aliased) «f3f*©1?, * ft 

5 ft . B««lBflT?*5J:5li:l*5, ^ o T , X y 7 1 2 T 18 SU £ ft ft B S ^ -fe 
/I/ © ffl * , X -r >y y 2 0 T 5 US BU £ ft ft , c©SfcoT^5>f^-b;KOfii:Jt8?ti, pi ^ 40 
© IB © 4> fc aft ft £ T © <K * -fe As # , X -r -y y l 2 1? m gij £ ft ft # * -tr ;U © ffi & 5? 9 fifc 
ft £ 0 dftfi, X r- >y y 2 2 T* ig C S o X f v ~f 2 4 T* , ^l?©A^6^^1«Jg^r, ^ft 

[ 0 0 2 1 ] 

0 5 ft , @4©7n-ft-K:^oTfffoni)Jl!!lCS!¥«$IT'$5o 
[ 0 0 2 2 ] 

0 6 ft, A^W^«^^#ffi^^MaM©iSftte<fcO : A^W=5:^Jfi^^ttL^^^S^© 
■ »**LT^5. & Hi m ft , f €H2 6, if 2 8 t5 t O'i 5 ^^fiSOta 3 0 
§A/"PV>§ 0 Stt a M © H © <f T* « , iy>^fi«©i«3 0!:t2 6i:©«S!!) , 1 Mm 
2 8 0tOtrai;«S«li:x^iJ7'>>y2hT^5 (alias) . fit, ^ftfCJE&CT 50 



(8) 



JP 2004-525715 A 2004. 8. 26 



<, fr-DlS?SSLTLf 7Ch&<> l& D & n 5 0 ± IE © ^ S ^ i* M fi . M « SE H 
[ 0 0 2 3 ] 

i-7 it. ±mi,rc&ffiic& zmmznfi? ztztbicm^z c tA'ftsso^iwnyif 

n - * 1 3 2 Ol H5%iLTV^5o nyifa-^ 1 3 2(i, f^ffllgll 3 4, 'J — K 

* >"J — ^ 1) 136, 7 V^ATi'-trX^t'J 1 38. A-Ff-fX^gll 4 0. 
^Xyi/^F7^^1 4 2 feJ;t>'f^ X^U-f 1 4 4, *-t:-K 1 4 Sfeilf^^X 1 10 

5 0%#-ri)3---tJ ? AtB^lE]SSl 4 6 * fll * T T , iT(i±tI©/5X 1 5 2^Mt8 
8ShT^§', 3 4(i, R O M 1 3 6 , RAMI 3 8 Sfttt^ - Kf-f X 
^ g B l 4 0rti:|g|ftJnT^§yn^7A^yx h7?->3 y^HfTLt, - £ ffi <d 
ffll^iRl, cn^rRAMl 38Sf:ttA-Ff^^lIl 4 0^tC^lft-r5Cli:^ 
T*#3„ ilBLfcififcT* — ^I«LT^T, to H fc: £> <^ T « . 0 4 tc ^ 
L , frOHUCfctfHStcR^WKSSSSftfcX^'y:/**??*-*. 7" a ^7* ^ A fi , tit * 

fi?«x.tfn>'>^hT f i'X^±{C^^LT, EfttiCttfTt*. £ fc ti , * y h "7 - * 

* 1 3 2 te: , jfiWfcn:/^— ^rn^AOlWISTTfi&ff-r**:*, J» * W K % ±I2L 20 
ft 8 ffi K i o T f-^^Mit^gf^^StS. ^gfi3y(?a-?132a, $ 

-^yn^7Aif B *ffl^TlffL, ±IE©<fc3fc«LSS*IH»-r*„ 
[0E©fS*;S:tSiW] 
[ 0 0 2 4 ] 

[0 i ] mte%%mmicfifc?&ffim*m?%mM*MffifticKLrcmT«&z><, 

[02] gilDI5i^'J7->y^©A^W4Ii^*lT^5|l cDM^MPtSWtc^ 
LftBTfeS. 

[0 3] ■«l«3T?a»t5ft«)OaSffiOB*iBOIK*«IIS»K^LfeHT?fe5. 
[04] B2fc*lftJ:7 4i«*'6AaW4ii4lR0»<8IIS**Lt^57n-^ 30 
-V - FT'$5, 

[05] 04OfilOft¥«a»T*S. 

[06] ±8BLfcAJftW*«3i©l»*O«C«T<0B«*^UfcH-e*S. 

[07] ± HE O ft « t «t 'S ffl S * * fT f S i& K ffl V S C t^T'f §I©^lW3>lfa 

[ «F # © IK B£ ] 
[ 0 0 2 5 ] 
2.26 #<D^$ 

4, 30 xe^^eioiRiii 

6 , 2 8 iflx 1 40 



(9) 



JP 2004-525715 A 2004. 8. 26 




I® 4 ] 



] 



| e pin s o s b ^.r^^tzJu^Dii ^aarga p ^ 12 
j d>8 Baxr^m^tz^tx^igDo gasrea 



D< B rdn V&&*?t2lKD&D L £Wft& 



-18 



IE 



-20 



24 




(10) 



JP 2004-525715 A 2004. 8. 26 



[0 7] 



134 , 136 / 138 



CPU 




ROM 




RAM 










t 



-132 



-152 



142 _L ^ 146 



I/O ^ 



144 



(11) 



JP 2004-525715 A 2004. 8. 26 



\^ > 7 \y V h ] 



(12) PiTERNAHONAI. APPLICATION WJBUSHKP UNDER THE PATENT COOPERATION TREATY fPCT) 



(19) World latelkctoaJ Property Organization 



[IDBIUlQIIQlSaBIDQillillQIIIl 



(43) Interna tloa*l Publication Date 
14 October 2001 (24,IOl1001) 



PCT 



WO 02/084594 A2 



(21) lateraaOoaal AppQcabeo dumber: KrTXHK2A»7CK 

(22) latarotianal raisg Date? 21 Rtstaoy 2002(21 

(25) nan l.inyi *a *: Uaglisb 



18 April 2001 (liOOWl) UB 

(71) A pplicaa t (/or eff .SkBej except U^. VOXAR 
LIMITED rOB/CTB]; Bcmdncton Ban]. 2 Amknoo 
Place, Edinburgh EH6 5NP (GB1. 

(72) Iww w n; *** 

. (75) lavcnrorsMppacaari Qcr USenfy): PAPAGCORCSOU, 



Parka [GKAUBj: 21/16 Blair Sow*. Hrftabcgh E1I11QR 
(OUX POOUC tea ItaWOBJ: 38 Prmce Oc^cm Street, lid- 
iaburxh L1H6 4AT ((H)). 

(74) Ag***: R4.INES. Miles Job*, d .1; D YUmg 4 Co, 21 
New Ferrer Lane, London E3C4A IDA (GB). 



(SI) Dc^bmhI States (mMionaQ: All AC, Al, AM, AT. AIJ, 
A7, RA. BB. Ba BR, BY. BZ.CA.CH.CN.OO.OI, CXI. 
CZ. DE, DK. DM. DZ. EC EE. ES, IT. OB. CD, GB. OH. 
GM. IK, 1IU. ID. IL. IN. IS. JP. EE, KG, KP. KB. EZ, LC. 
LE.LK.LS.Lr. LU. LV. MA, MD, MU, MK. MN. MW. 
MX. MX, NO, NV, OM. PH. PI, PI*. IK). KU. SI3. SB, SO. 
SL SK, SI, TJ. TM. TN. TR, TT. T7, I'A. IJG. OS, U7, 
VS. YU. ZA, ZM. ZW. 



(84) Designated States frtjrioaaOj AR1PO pttcm (Gil, GM. 
Kl!, 1-S, MW. MZ. SO, SI, S7, T7, IJG. ZM. ZW). 
Ibundan patent (AM, A7, BY, KG. KZ, Ml), K1I, TJ. TM X 



{ (54) TMK CORRBCTION OP BOUNDARY ARTEFACTS IN IMAGE DATA PROCESSING 



(57) Abstract: 



3 



gjs I Capture voxels characterised by density D 1 
— | Identify set of voxels LV with 3™ £D ^Bma, 



Identify set of voxels Ou with D > B™u 



Identify set of voxels L\ with D < ~\j 



in 



Apply approximatd sphericaJ morphological 
dilation to (X and Ol to form D'„ v D' L 



Identify set Dp, given by intersection 
of r> 0 and D'l 



JL 



-20 



Subtract Opv from Due to form Out 
; 

Display D*ul 



4 



3^- 



in which cats of im»ys elrmenU 
having display values ouukic of 
■ (argot renfo B of display values 
tic c iL ' n respectively nwpfn>lfl(g» 
►caTTy dilated. The i ataa a cti ow 
bciuccn the morphologically 
dilaicij sets of image dcxPCO ti 



ftum the «el of lm 
having the target range B of 
display rabies. This remotes 
image elements taccrrccdy 
appearing to have ditpfcry nines 
corresponding <*> the largel 



having 

display values either tide of 
the target range B of imaging 
may he two-dimensional or 
thtcc-dimciuiooal imaging. 
Tho morphologically dilatation 
is preferably performed 



al 

ml hiiving a radius of 



i 



(12) 



JP 2004-525715 A 2004. 8. 26 



WO 02/084594 A2 IBBIMMjljligMjgfflMM 



Utropcrs pens CAE Bli. CH. CX DC DK_ US. I T. I-H, t~ur t^-U^r crxiu oO*, iabsitx. rtfir u X^U- 

CiB. CIK, U', rr. Ul, MC, NU FT, SI', TR\ QAM p«eu an AVtaf rwCoAjcnrf^jxSrroaioni' ^iw^u c « «** 

fRP. BJ. CF, CO, d. Of. OA, ON. GO, GW. Ml, MR. tmj efcachrcgtlr torn tfthe PCT Germ*. 
SB. SN, TD. TG). 

PHbfidMdb 

tritkemi BamKOodJ i«a * cA (xpo-f corf to 6* rrpabiis'trd 
upon rtetipt of that report 



(13) JP 2004-525715 A 2004.8.26 



WO 02JOUS94 PCT/GB01/00768 
TITLE OP THE INVENTION 

CORRECTION op BOVNPARY ARTEFACTS^ 
IMAGE DATA PROCESSING 

5 

FIELD OP THE INVENTION 

This invention relates to the field of image data processing. More particularly, 
this invention relates to the display of image elements having display values lying 
10 within a range of display values located within a broader overall range of display 
values. 

BACKGROUND ART 

15 It is known to provide image display systems in which image elements, 

whether they are representing a two-dimensional image or a three-dimensional image, 
are displayed to a user with a selection and modified processing of at Least some of the 
image elements being made in dependence upon the display values for those image 
elements. In certain cases, these image elements are directly related to a physical 

20 characteristic of an object through the' measurement of signal strength values from a 
defection device. As an example, in medical imaging applications signal values 
representing the signals returned from CAT scanning or MRI scanning may be 
displayed with display elements having their intensity or colour controlled by the value 
of the signal returned for each particular image element In order to improve the ease 

25 of interpretation of such images it is known to map different colours to different ranges 
of display value such that particular features, e.g. blood vessels, may be made more 
visible within the image. 

A particular problem arises when it is desired to selectively process or display 
those image elements falling within a range of display values that is itself located 

30 within a broader full range of display values for the image as a whole. Figure 1 of the 
accompanying drawings illustrates such an image. The image contains a region 2 of 
image elements having high display values, such as a region of bone within a medical 
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image. This region 2 having high display values is located within an overall 
background 4 of image elements having low display values, such as image ele men ts 
r^ pr ^ y r^pg soft tissue. Also present within this image is a blood vessel 6 containing 
a contrast enhancing ftfc^ 1 it thai is represented by image elements having display values 
5 Located somewhere between the high values of the bone region 2 and the low values of 
the background soft tissue region 4. 

Figure 2 of the accompanying drawings illustrates the problem that occurs due 
to the finite resolution of real imaging systems in that at the interface between the bone 
region 2 and the soft tissue 4 there is a boundary formed of pixels having values tying 

10 somewhere between the high values of the bone region 2 and the low values of the soft 
tissue region 4. In many cases, these boundary pixels will have display values closely 
similar to the display values corresponding to the blood vessel 6. Accordingly, if the 
image of Figure 2 is processed to enhance the appearance, or otherwise select in some 
way, the display elements having display values corresponding to the blood vessels 6, 

15 then this will erroneously also highlight a region at the interface between the bone 
region 2 and the soft tissue 4 that due to aliasing effects appears to have the 
appropriate display value. 
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SUMMARY OF THE INVENTION 

Viewed from one aspect the present invention provides a method of processing 
image data formed of an array of image elements, each image dement having a display 
value, said method comprising the steps of: 
5 identifying a set of target image elements having display values withra a target 

range of display values from B^*, to B^; 

identifying a set of low display value image elements having display values 
below T^; 

identifying a set of high display value image elements having display values 
10 above Tnaxj 

applying a morphological dilation to said set of low display value image 
elements to generate a dilated set of low display value image elements; 

applying a morphological dilation to said set of high display value image 
elements to generate a dilated set of high display value image elements; 
15 identifying an intersection set of image elements that are present in both said 

dilated set of low display value image elements and said dilated set of high display 
value image elements; and 

removing from said set of target image elements any image elements also 
present within said intersection set of image elements to form a modified set of target 
20 image elements. 

The invention recognises that the image elements that are incorrectly appearing 
to have display values corresponding to the target range of display values occur at the 
interfaces between regions of image elements with display values at opposite sides of 

25 the target range and if these interface image elements can be reliably identified, then 
they may be removed. The invention identifies such interface image elements by 
morphologically dilating the regions (i.e. sets of image elements) at either side of the 
target region and then determining the intersection of those dilated regions. The 
dilated regions will intersect where they touch one another and accordingly this 

30 technique selectively identifies the interface regions rather than picking up image 
elements that lie at an interface between the target region and a region of image 
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dements having display values outside of the target range. When these intrrfarr 

elements having a display value within the target range and the modified set of image 
elements having the target display value can then be displayed with reduced interface 
S aliasing artefacts ox in some other way processed, e.g. volume measured. 

The technique of the invention is able to identify and remove the image 
elements producing the erroneous artefacts within the image without significantly 
impacting the display of image elements that are correctly iden t ified. By "comparison, 
a technique that merely morphologically expanded either the high display value region 
10 or the low display value region to overwrite any interface image elements would also 
be likely to overwrite the edges or the fine detail in the display of the regions image 
elements that correctly have the target range of display values, e.g. fine capillaries may 
be lost from the image. 

It will be appreciated thai the technique described above could be applied to 
15 both two dimensional images and three dimensional images. The invention is 
particularly well suited to use within imaging arrays of three dimensional voxel data as 
such images can be significantly degraded by the aliasing artefacts between regions as 
discussed above. 

The morphological dilatation could take a variety of forms depending upon the 
20 particular circumstances, but preferably has the form of a spherical morphological 
dilatation whereby each voxel is projected onto all the points within a quasi spherical 
surrounding region. 

The spherical structuring element used in this morphological dilatation could 
have a variety of sizes, but the invention has been found to be particularly effective 
25 when the structuring element has a radius of between 2 and 3 voxel sizes, and more 
preferably substantially 2.5 voxel sizes. 
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that correspond to the ranges being identified and are used to control the way in which 
those image elements are displayed. In preferred practical systems, the display values 
do not correspond to what would normally be regarded as visual properties, such as 
5 colour or intensity, but instead relate to detected signal values from measuring systems 
such as CAT scanners, and MRI scanners, ultrasound scanners end PET systems. 

The invention is particularly well suited to the removal of artefacts from 
images when attempting to discriminate between blood vessels containing contrast 
enhancing agents, soft tissue and bone within a medical diagnostic, such as angiogram, 
10 image. 

Other aspects of this invention provide apparatus for processing image data and 
a computer program for controlling a computer to process image data in accordance 
with the above described techniques. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

An embodiment of the invention will now be described, by way of example 
only, with reference to the ecconrpenying drawings in which: 

Figure 1 schematically illustrates an image having regions corresponding to 
5 different display values; 

Figure 2 schematically illustrates the image of Figure 1 showing the aliasing 
artefacts *h m arise; 

Figure 3 schematically illustrates the selection of a target range of display 
values to highlight within an image; 

10 Figure 4 is a flow diagram illustrating a technique for removing artefacts from 

an image such as that shown in Figure 2; 

Figure 5 is a mathematical representation of the processing of Figure 4; 

Figure 6 illustrates an image bom before and after the removal of the artefacts 
described above; and 

15 Figure 7 schematically illustrates a general purpose computer of the type that 

may be used to carry out processing in accordance with the above techniques. 
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DETAILED DESCRIPTION 

Image da t a, such as a collection of say 100 twcMiimerisional 512*512 images 
copBCted from a CAT scanner, a MR1 scanner, an ultrasound scanner or a PET 
(Positron emission tomography system), may be subject to image proces sin g in 
5 accordance with known techniques to" produce a thxee-dimensianal representation of 
the structure imaged (various user selected two-dimension projections of the three- 
dimensional representation are typically displayed an a computer monitor). The 
techniques for £cxxci sting such 0 u* f^r^lirn cnston^l rcpr cs cntttio tts of structures from 
collections of two-dimensional images are known in themselves and will not be 
10 described further hereirL 

In order to increase the uoderstandabOity of the threc^iimensioTial 
representations generated it is known to select ranges of display values for highlighting 
or selective processing in some other way. As an example, voxels having a particular 
range of display value may be tinted with a vivid colour to stand out within the image 
15 or may be selected for removal from the image to reveal other more interesting 
features. 

Figure 3 is a histogram illustrating the frequency of occurrence of voxels 
within an image as a function of the signal value associated with those voxels. As 
illustrated in Figure 3, there is a target range B of interest that may have its upper B^m 

20 and lower Bmb boundaries selected by the user The user may select this range to try 
to pick out blood vessels containing a contrast enhancing agent, such as maybe desired 
when performing angiography. This target region is positioned within the broader full 
range of signal values that may be generated with a lower signal value range A 
bounding the target range B on its lower side. Similarly, a higher value range C 

25 bounds the target range B on its higher side. As previously discussed, a problem can 
arise when tissue that returns a signal value within the range A abuts tissue that returns 
a signal value within the range C as the spatial resolution of the system and the 
consequences of finite resolution sampling may result in an image element, voxel, 
being generated at this interface with a signal value corresponding to the range B even 

30 though this is not truly a blood vessel region that it is desired to highlight. 
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Figure 4 Is a flow diagram fflustrating the artefact removal technique. At step 
10, the voxels that are to form the image are captured with each one having a 
corresponding display value, in this case the display values are representative of signal 
5 strength values from a detection device, and indicate a density value D. At step 12, the 
set of voxels having a density within the user specified target range B of density values 
illustrated in Figure 3 is identified. At step 14, the set of voxels having a density 
corresponding to the range C of Figure 3 is identified. At step 16, the set of voxels 
having a density corresponding to the range A of Figure 3 is identified. 

10 The sets of voxels identified at steps 14 and 16 are then subject to respective 

spherical morphological dilatation to produce dilated sets of voxels. The 
morphological dilatation can take a variety of forms, but in this example uses a quasi- 
spherical structuring element based upon a voxel approximation to a sphere having a 
radius of between 2 and 3 voxels, but preferably substantially 2.5 voxels. The 

IS morphological dilatation with such a spherical structure acts to project each voxel 
value onto all of the voxels within a region surrounding the starting voxel defined by 
the spherical structure. This slightly expands/dilates the region concerned. 

At step 20, the two dilated sets of image elements are compared to identify 
image elements appearing within both sets. These image elements correspond to the 

20 interface regions between the two sets. These interface regions are the place where 
voxels incorrectly aliased so that they appear to be within the target range B may 
appear. Accordingly, the set of target voxels identified at step 12 is compared with this 
intersecting set of voxels identified at step 20 and any voxels that appear in both sets 
are removed from the set of voxels identified at step 12. This occurs at step 22. At 

25 step 24, the resulting set of target voxels having the interface artefact removed from 
them are displayed. 

Figure 5 is a mathematical representation of the processing that occurs in 
accordance with the flow diagram of Figure 4. 
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Figure 6 illustrates o Before image in which the artefact is present end en After 
image in which the artefact is not present. Each image contains a region of bane 26, 
blood vessels 28 end a soft tissue region 30. In the Before image, the interface 
between the soft tissue region 30 and the bone 26 aliases to display values similar to 
5 those of the blood vessels 28 and accordingly this interface region is tinted in the same 
way as the blood vessel 28 making the image more difficult to interpret. The After 
image shows the effect of applying the above technique to removing this artefact. This 
artefact is removed without inappropriately changing the shape of the regions 
appearing within the image and without the possibility of obliterating fine detail that 
1 0 may occur by simply dilating one region outside of the target range. 

Figure 7 schematically illustrates a general purpose computer 132 of the type 
that may be used to perform processing in accordance with the above described 
techniques. The computer 132 includes a central processing unit 134. a read only 
memory 136, a random access memory 138, a hard disk drive 140, a display driver 142 

15 and display 144 and a user input/output circuit 146 with a keyboard 148 and mouse 
ISO all connected via a common bus 152. The central processing unit 134 may 
execute program instructions stored within the ROM 136, the RAM 138 or the hard 
disk drive 140 to carry out processing of data values that may be stored within the 
RAM 138 or the hard disk drive 140. Data values may represent the image data 

20 described above and the processing may carry out the steps illustrated in Figure 4 and 
as expressed mathematically in Figure 5. The program may be written in a wide 
variety of different programming languages. The computer program itself may be 
stored and distributed on a recording medium, such as a compact disc, or may be 
downloaded over a network link (not illustrated). The general purpose computer 132 

25 when operating under control of an appropriate computer program effectively forms an 
apparatus for processing image data in accordance with the above described technique. 
The general purpose computer 132 also performs the method as described above and 
operates using a computer program product having appropriate code portions (logic) 
for controlling the processing as described above. 
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CLAIMS 

I • A meth od of processing image data formed of an array of image elements, each 
image element having a display value, said method comprising the steps of: 
5 identifying a set of target image elements having display values within a target 

range of display values from Bgfa to 

identifying a set of low display value image elements having display values 

below B«jja; 

identifying a set of high display value image elements having display values 
10 above Bob 

applying a morphological dilation to said set of low display value image 
elements to generate a H»*tw< set of low display value image elements; 

applying a morphological dilation to said set of high display value image 
elements to generate a dilated set of high display value image elements; 
15 identifying an intersection set of image elements that are present in both said 

dilated set of low display value image elements and said dilated set of high display 
value image elements; and 

removing from said set of target image elements any image elements also 
present within said intersection set of image elements to form a modified set of target 
20 image elements. 

2. A method as claimed in claim 1, wherein said image elements are voxels and 
- said array is a three dimensional array of voxels. 

25 3. A method as claimed in claim 2, wherein said morphological dilation applied to 
said set of low display value image elements uses a quasi-spherical structuring 
element. 

4. A method as claimed in any one of claims 2 and 3, wherein said morphological 
30 dilation applied to said set of high display value image elements uses a quasi-spherical 
structuring element- 



to 
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5. A method as claimed in any one of claims 3 and 4, wherein said quasi-spherical 
structuring element has a radius of between 2 and 3 voxels. 

6. A method as claimed in claim 5, w heron said quasi •spherical structuring 
5 element has ft radius of substantially 2.5 voxels. 

7. A method as claimed in any one of the preceding claims, wherein said display 
values represent a measured signal strength value returned from a portion of a subject 
being imaged mapped to a corresponding image element. 

10 

8. A method as claimed in claim 7, wherein said measured signal strength values 
arc detected using one of: 

CAT scanning; and 
MR1 scanning; 
IS ultrasound scanning and 

PET. 

9. A method as claimed in any one of the preceding claims, wherein said set of 
low value image elements correspond to soft tissue, said set of high value image 

20 elements correspond to bone and said target set of image elements correspond to blood 
vessels containing a contrast enhancing agent 

10. Apparatus for processing image data formed of an array of image elements, 
each image element having a display value, said apparatus comprising: 

25 a target set identifier operable to identify a set of target image elements having 

display values within a target range of display values from B min to Bm„; 

a low display value set identifier operable to identify a set of low display value 
image elements having display values below B^; 

a high display value set identifier operable to identify a set of high display 
30 value image elements having display values above Bm,; 
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a low set dilator operable to apply a morphological dilation to said set of low 
display value image elements to generate a dilated set of low display value image 

a high set dilator operable to apply a morphological dilation to said set of high 
5 display value image elements to generate a dilated set of high display value image 
el eM^n t^ 

an intersection identifier operable to identifying an intersection set of image 
elements that are present in both said dilated set of low display value image elements 
and said dilated set of high display value image dements; and 
10 an intersection remover operable to remove from said set of target image 

elements any image dements also present within said intersection set of image 
el e m ents to form a modified set of target image dements. 

11. Apparatus as claimed in claim 10, wherein said image elements are voxels and 
1 5 said array is a three dimensional array of voxels. 

12. Apparatus as claimed in claim 1 1, wherein said morphological dilation applied 
to said set of low display value image dements uses a quasi-spherical structuring 
element. 

20 

13. Apparatus as claimed in any one of claims 11 and 12. wherein said 
morphological dilation applied to said set of high display value image elements uses a 
qua si -spherical structuring element 

25 14. Apparatus as claimed in any one of claims 12 and 13, wherein said quasi- 
spherical structuring element has a radius of between 2 and 3 voxels. 

15. Apparatus as claimed in claim 14. wherein said quasi-spherical structuring 
dement has a radius of substantially 2.5 voxels. 

30 
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16. Apparatus as claimed in any one of claims 10 to 15, where in said display 
values represent a measured signal strength value returned from a portion of a subject 
being imaged mapped to a corresponding image element 

5 17. Apparatus as claimed in claim 16, wherein said measured signal strength 
values arc detected using one of: 
CAT seaming; 
MRI scanning; 

10 PET. 

18. Apparatus as claimed in any one of claims 9 to 17, wherein said set of low 
value image elements correspond to soft tissue, said set of high value image elements 
correspond to bone and said target set of image elements correspond to blood vessels 

IS containing a contrast enhancing agent . 

19. A computer program product bearing a computer program operable to control a 
computer to process image data formed of an array of image elements, each image 
element having a display value, said computer program comprising: 

20 target set identifying logic operable to identify a set of target image elements 

having display values within a target range of display values from to IW; 

low display value set identifying logic operable to identify a set of low display 
vahic image elements having display values below B m jr ; 

high display value set identifying logic operable to identify a set of high 
25 display value image elements having display values above B^; 

low set dilating logic operable to apply a morphological dilation to said set of 
low display value image elements to generate a dilated set of low display value image 
elements; 

high set dilating logic operable to apply a morphological dilation to said set of 
30 high display value image elements to generate a dilated set of high display value image 
elements; 
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interaction irimti tying logic operable to identifying an intersection set of 
image elements that are present in both said dilated set of low display value image 
elements and said dilated set of high display value image elements; and 

5 elements any image elements also present within said intersection set of image 
elements to form a modified set of target image elements. 

20. A computer program product as claimed in claim 19, wherein said image 
elements are voxels and said array is a three dimensional array of voxels. 

10 

21. A computer program product as claimed in claim 20, wherein said 
morphological dilation applied to said set of tow display value image dements uses a 
quasi-spherical structuring element 

IS 22. A computer program product as claimed in any one of claims 20 and 21, 
wherein said morphological dilation applied to said set of high display value image 
elements uses a quasi -spherical structuring element 

23. A computer program product as claimed in any one of claims 21 and 22, 
20 wherein said quasi-spherical structuring element has a radius of between 2 and 3 

voxels. 

24. A computer program product as claimed in claim 23, wherein said quasi- 
spherical structuring element has a radius of substantially 2.5 voxels. 

25 

25. A computer program product as claimed in any one of claims 19 to 24, wherein 
said display values represent a measured signal strength value returned from a portion 
of a subject being imaged mapped to a corresponding image element. 

30 26. A computer program product as claimed in claim 25, wherein said measured 
signal strength values are detected using one of: 
CAT scanning; 
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